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ABSTRACT 

Coupling  the  3-d  ocean  optical  imagery  with  3-d  circulation  models  provides  a  new  capability  to  understand  coastal 
processes.  Particle  distribution  derived  from  ocean  color  optical  properties  were  coupled  with  numerical  circulation 
models  to  determine  a  24  hour  forecast  of  particle  concentrations. 

A  3-d  particle  concentration  field  for  the  coastal  ocean  was  created  by  extending  the  surface  satellite  bio-optical 
properties  vertically  by  parameterzing  an  expediential  Gaussian  depth  profile.  The  shape  of  the  vertical  particle  profile 
was  constrained  by  1)  the  depth  of  the  1%  light  level  2)  the  mixed  layer  depth  3)  the  intensity  of  the  layer  stratificatior 
4)  and  subsurface  current  field  and  the  surface  bio-optical  properties.  These  properties  were  obtained  from  MODIS 
ocean  optical  products  (phytoplankton  absorption  and  backscattering)  and  the  Intra-America  Sea  Nowcast  Forecast 
System  -  Naval  Coastal  Ocean  Model. 

The  3-d  particle  distribution  was  imbedded  into  a  3-d  circulation  model  and  the  particles  advected  hourly  using  foreca 
model  3-d  current.  The  particles  were  diffused,  dispersed  and  differentially  settled  during  the  advection  processes. 
Following  the  24  hour  advection,  the  resultant  particle  distribution  were  accumulated  into  1  km  spatial  grid  and 
vertically  to  a  1  attenuation  length  (satellite  penetration  depth)  and  the  forecast  ocean  color  backscattering  image 
determined.  The  forecast  image  was  compared  with  the  next  day  ocean  color  backscattering  image  to  define  the  error 
budget. 

The  ocean  color  particle  tracking,  defines  fine  spatial  scales  processes  such  as  local  upwelling  and  downwelling,  whic 
are  essential  in  understanding  the  coupling  of  physical  and  bio-optical  processes.  The  methods  provide  new  capability 
for  characterizing  how  subsurface  particles  layers  change  in  response  to  cross  and  along  shelf  exchange  processes. 
Results  show  methods  to  forecast  satellite  optical  properties  in  coastal  areas  and  examine  how  sequential  MODIS 
imagery  of  the  particle  scattering  is  related  to  particle  transport  and  physical  processes 
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1.  INTRODUCTION 

Remote  sensing  of  ocean  color  is  used  to  characterize  and  monitor  coastal  processes.  Our  capability  to  exploit  ocean 
color  to  retrieve  the  near  surface  bio-optical  properties  is  limited  to  within  the  upper  first  attenuation  length.  The  vertical 
variability  of  the  bio-optical  properties  changes  significantly  in  response  to  both  biological  and  physical  processes. 

These  subsurface  biological  processes  are  important  to  understand  the  optical  layers  such  as  the  location  and  intensity  of 
the  chlorophyll  maximum  and  bottom  resuspension  layers.  These  layers  influence  the  shape  of  the  vertical  profile  which 
is  tied  to  the  near  surface  observation  of  surface  satellite  imagery1. 


Coastal  Ocean  Remote  Sensing,  edited  by  Robert  J.  Frouin,  ZhongPing  Lee, 
Proc.  of  SPIE  Vol.  6680,  66800S,  (2007)  0277-786X/07/$18  •  doi:  1 0. 1 1 1 7/1 2.737201 


Proc.  of  SPIE  Vol.  6680  66800S-1 


20080211225 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 

gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  the  Department  of  Defense,  Executive  Services  and  Communications  Directorate  (0704-0188)  Respondents  should  be  aware 
that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB 
control  number 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

0 1  -02-2008  Conference  Proceeding 

3.  DATES  COVERED  (From  -  To ) 

4.  TITLE  AND  SUBTITLE 

Forecasting  Coastal  Optical  Properties  Using  Ocean  Color  and  Coastal 
Circulation  Models 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

0601 153N 

6.  AUTHOR(S) 

Robert  A.  Arnone,  Brandon  J.  Casey,  Dong  Shan  Ko,  Peter  M.  Flynn,  L. 
Carrolor,  Sherwin  Ladner 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

73-8404-B6-5 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

Oceanography  Division 

Stennis  Space  Center,  MS  39529-5004 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NRL/PP/7330-07-724 1 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 

800  N.  Quincy  St. 

Arlington,  VA  22217-5660 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

ONR 

11.  SPONSOR/MONITOR'S  REPORT 

NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Coupling  the  3-D  ocean  optical  imagery  with  3-D  circulation  models  provides  a  new  capability  to  understand  coastal  processes. 
Particle  distribution  derived  from  ocean  color  optical  properties  were  coupled  with  numerical  circulation  models  to  determine  a  24 
hour  forecast  of  particle  concentrations.  A  3-D  particle  concentration  field  for  the  coastal  ocean  was  created  by  extending  the 
surface  satellite  bio-optical  properties  vertically  by  parameterzing  an  expediential  Gaussian  depth  profile.  The  shape  of  the  vertical 
particle  profile  was  constrained  by  1)  the  depth  of  the  1%  light  level  2)  the  mixed  layer  depth  3)  the  intensity  of  the  layer 
stratification  and  4)  subsurface  current  field  and  the  surface  bio-optical  properties  were  obtained  by  MODAS  ocean  optical  products 
(phytoplankton  absorption  and  backscattering)  and  the  Intra-Americal  Sea  Nowcast  Forecast  System-  Naval  Coastal  Ocean  Model.... 

15.  SUBJECT  TERMS 

ocean  optics,  satellites,  ocean  color,  ocean  models 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

ABSTRACT 

18.  NUMBER 

OF 

PAGES 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

Unclassified 

Unclassified 

Unclassified 

UL 

8 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Robert  A.  Arnone 


19b.  TELEPHONE  NUMBER  ( Include  area  code) 

228-688-5268 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std  Z39  18 


New  methods  are  required  to  extend  the  near  surface  remote  sensing  to  characterize  the  subsurface  properties.  These 
capabilities  are  important  for  understanding  the  total  exchange  processes  across  front  zones  and  the  transport  studies 
from  coastal  to  open  ocean  waters. 

Remote  sensing  ocean  color  provides  a  unique  capability  to  monitor  the  bio-optical  conditions  within  the  coastal  zone  by 
providing  a  ‘Nowcast”  of  the  current  conditions.  However,  coastal  managers  require  a  capability  to  “forecast”  the  bio- 
optical  conditions  to  be  used  for  decision  making.  Improved  methods  are  required  to  forecast  the  24  hour  surface  bio- 
optical  properties  and  can  be  accomplished  by  combining  remote  sensing  products  with  forecast  circulation  models  2  \ 

The  near  surface  particle  distribution  field  can  be  retrieved  from  the  remote  sensed  backscattering  intensity  4,5 .  This 
particle  field  is  coupled  with  modeled  currents  to  advect  the  particle  field  to  provide  a  24  hour  forecast.  The  particle  field 
distribution  from  the  satellite  backscattering  imagery  provide  the  initial  “seed”  field  for  the  model  currents  as  a 
representation  of  the  spatial  concentrations  and  distribution  of  particles.  The  model  currents,  transports  these  particle 
trajectory  on  hourly  time  steps  to  determine  the  24  particle  concentration.  The  particles  are  than  accumulated  into  1  km 
bins  representing  similar  resolution  as  ocean  color  satellites  and  converted  back  to  a  backscattering  coefficient.  This 
particle  tracking  assumes  particle  settling  and  diffusion  and  only  particle  “seeding”  at  the  surface.  (The  previous  3-d 
vertical  bio-optical  model  will  be  used  to  generate  a  3-d  particle  “seed”  for  the  circulation  in  the  near  future.)  Biological 
processes  such  as  plankton  growth  and  decay  as  not  included  in  this  approach.  This  is  an  initial  effort  to  estimate  the 
particle  forecast  to  assess  how  much  of  the  variability  can  be  assessed  based  on  simply  advection  processes.  The 
forecast  capability  of  the  particle  field  provides  a  useful  tool  to  assess  the  changing  coastal  conditions  and  provide 
decision  making  for  the  future. 


2.  DETERMINING  THE  VERTICAL  3-D  OPTICAL  VARIABILITY 


The  extensions  of  the  satellite  optical  properties  were  performed  by  characterizing  the  shape  of  the  vertical  as  a 
Gaussian1 ,2.  The  vertical  profiles  of  bio-optical  profiles  such  as  the  chlorophyll  maximum  layer  have  been  shown  to  be 
well  defined  using  a  Gaussian  shape  6.  We  developed  new  methods  to  provide  constraints  on  the  profile  shape  which  are 
controlled  by  the  link  between  the  mixed  layer  depth  (MLD),  intensity  of  the  mixed  layer  (IMLD)  and  the  1  percent  light 
level.  These  constrains  were  determined  by  analyses  of  many  profiles  of  temperature,  salinity,  fluorescence  and 
backscattering  in  a  wide  variety  of  open  and  coastal  waters  representing  a  variety  of  water  types  and  process.  In  this 
study,  the  physical  properties  of  MLD  and  IMLD  are  provided  through  the  physical  circulation  model  and  the  subsurface 
light  fields  of  the  1  %  light  level  are  estimated  using  the  surface  satellite  inherent  optical  properties7 .  These  shape 
parameters  of  the  Gaussian  are  than  linked  to  the  surface  satellite  surface  chlorophyll  1  .  The  vertical  Gaussian  shape  is 
defined  by: 
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where  parameters  are  defined  in  Figure  1  and  (k)  is  the  decay  of  the  background  chlorophyll.  Adjustment  of  these 
parameters  allows  various  vertical  shapes  to  be  defined.  The  shape  of  the  profiles  is  coupled  to  the  surface  chlorophyll 
by  setting  the  integration  of  the  1st  optical  depth  of  this  shape  to  the  satellite  chlorophyll  value. 
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Figure  1  -  A.  Shape  of  the  Bio-optical  profile  identifying  the  components  (Zm,a,  h,  Bo).  B.  The  location  of  the 
optical  layer  constrained  between  the  mixed  layer  depth  from  the  density  profile  and  the  1%  light  level.  The 
intensity  of  density  stratification  defines  the  height  or  elevation  of  the  Gaussian  peak. 


The  parameters  which  are  used  to  constraint  this  profile  shape  are  obtained  from  both  the  surface  satellite  optics  and  the 
Navy  Coastal  Ocean  Model  (NCOM)8.  The  density  field  from  NCOM  is  used  to  determine  the  MLD  and  the  MLD. 
NCOM  is  a  41  layer  model  currently  used  for  Navy  operations.  The  models  is  forced  with  COAMPS  winds  and  obtains 
boundary  conditions  fro  the  “Global  NCOM”.  These  models  assimilate  the  sea  surface  height  fields  from  satellite 
altimeter  and  sea  surface  temperature  from  AVHRR  9.  The  vertical  density  field  from  NCOM  is  used  to  estimate  the 
MLD  and  the  density  derivative  is  used  to  estimate  the  degree  of  stratification  or  the  IMLD.  Figure  lb  illustrates  how  the 
location  of  the  optical  layer  is  linked  to  the  MLD  and  the  IMLD.  NCOM  examples  of  the  MLD  and  the  IMLD  are 
shown  below. 


Figure  2.  Output  from  the  NCOM  vertical  density  profile  is  used  to  determine  the  Mixed  layer  depth  and  the 
intensity  of  the  mixed  layer  depth.  These  parameters  are  linked  with  the  vertical  bio-optical  properties.  Example 
of  the  NCOM  for  the  Mississippi  River  discharge  is  shown. 
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The  IMLD  indicate  the  degree  of  stratification  and  this  controls  the  sharpness  of  the  chlorophyll  layer  (a,  h).  The  1 
percent  light  level  depth  (PD)  was  computed  using  the  satellite  retrieved  absorption  and  backscattering  properties  as 
described  by7 10. 

The  depth  of  the  chlorophyll  max  (Zm)  is  constrained  between  the  MLD  and  the  PD  using  the  equation: 

Zm  =  f  (  MLD  +  (PD  -  MLD)  *  e  ('C1*IMLD)  )  (2) 

Similarly  the  Chlorophyll  Maximum  for  the  profiles  is  estimated  by: 

Chl-Max  =  f(  MLD  /  PD  *  C2/SfcChl)  (3) 

The  Gaussian  Spread  is  estimated  using: 


ct  =  f(  (C3-1)*EXP(-C4*IMLD  )  +  l)  (4) 


The  coefficients  Cl,  C2,  C3,  and  C4  were  determined  through  optimization  based  on  insitu  data  stations  of  the  vertical 
density  profiles  and  the  chlorophyll  profiles.  The  spatial  and  temporal  variability  of  these  coefficients  is  an  active  area  of 
research  which  includes  characterizing  the  response  of  these  coefficients  to  different  processes  which  include  bio-optical 
and  physical  coupling. 

The  results  of  the  3-d  model  are  shown  by  comparison  with  insitu  observations  of  fluorescence  and  density. 
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Figure  3:  A.  The  results  of  the  3-d  model  are  shown  in  comparison  with  insitu  density  and  the  fluorescence 
profile.  The  profile  shape  is  coupled  with  the  surface  satellite  chlorophyll  to  the  first  optical  depth.  The 
coefficients  for  this  profile  are  shown  for  Cl,  C2,  C3,  and  C4. 


B.  3-d  chlorophyll  volume  of  the  Northern  Gulf  of  Mexico.  The  surface  MODIS  chlorophyll  is  extended 
vertically  using  the  3-d  model  by  combining  with  the  NCOM  physical  properties.  Mississippi  River  discharge  is 
on  the  left.  Note  layers. 
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3.  FORECASTING  THE  OPTICAL  ENVIRONMENT 


The  coupling  of  remote  with  numerical  forecast  models  permits  the  capability  to  forecast  the  ocean  optical  environment. 
The  ocean  color  imagery  from  MODIS  has  been  used  to  determine  the  near  surface  backscattering  coefficient  (bb555). 
The  bb555  is  directly  related  the  number  of  particles  in  the  surface  waters  in  addition  to  the  particle  size,  composition 
and  index  of  refraction  5,11 .  We  assume  a  linear  relationship  between  the  particle  concentration  and  the  bb555  to  estimate 
the  surface  particle  numbers  in  a  satellite  image.  This  surface  particle  field  is  characterized  as  the  “seed”  field  of  the 
particle’s  spatial  distribution  of  surface  waters.  This  “seed”  field  is  inserted  into  the  surface  circulation  field  of  the 
NCOM  model  and  the  particles  are  dispersed  and  settled  at  hourly  time  sets  based  on  the  model  hour  forecast  field  for  a 
period  of  24  hours.  Note  that  the  satellite  particle  field  represents  only  near  surface  layer  and  not  the  3-d  particle  field  (in 
its  present  configuration.) 

Following  the  25  hour  particle  dispersion,  the  particles  are  accumulated  into  1  km  regional  pixels  and  transformed  back 
into  the  bb550.  This  resulting  25  hour  forecast  represents  the  forecast  backscattering  MODIS  image  which  can  be 
compared  with  the  next  day  MODIS  bb550  image.  The  difference  fields  represent  the  errors  associated  with  1)  model 
errors  in  circulation  2)  errors  is  particle  dispersion  3)  non-conservative  particle  interaction  including  settling  and  grown, 
accumulation  etc.  4)  particle  upwelling  and  downwelling. 
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Figure  4.  The  MODIS  backscattering  product  is  used  to  “Seed”  the  surface  model.  The  particles  are  forecast  24 
and  48  hourly  and  than  accumulated  into  a  forecast  bb555  image.  The  forecast  are  compared  with  the  next  day 
MODIS  bb555  imagery  to  confirm  errors. 


Present  issues  with  the  forecast  capability  are  linked  to:  1)  the  initial  seed  field  2)  inaccuracies  with  the  coastal  inputs 
(rivers  etc.)  3)  surface  convergence  and  divergence  fields.  The  errors  in  the  initial  seed  field  are  associated  with  non- 
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complete  satellite  coverage  from  clouds  etc.  These  non-complete  particle  coverage  results  in  particle  void  regions  which 
are  advected  in  the  forecast.  These  “voids”  regions  are  presently  being  filled  using  a  triangular  bb55  interpolation 
procedure12  which  provides  a  method  to  initialize  with  complete  coverage. 

Similar  problems  with  particle  void  areas  occur  when  as  particles  are  removed  with  no  replacement.  These  areas  are 
evident  at  river  months  where  the  initialized  particle  field  is  advected  away  from  the  source  leaving  a  void  (see  Figure  4 
-Arrow  pointing  to  River  plume).  As  these  particle  void  areas  is  generated,  they  are  observed  emanating  along  coastal 
areas.  These  are  can  /  should  be  filed  as  either  point  source  particle  field  (i.e.  River  particle  flux)  or  upwelled  particles. 
(This  river  flux  was  filled  in  Figure  5.) 

An  example  of  the  Mississippi  River  plume  disperion  on  to  the  Northern  Gulf  of  Mexico  is  shown  for  April  23,2007  and 
forecasting  April  24,  2007. 


0  0005  0  002  0.0079  0  032  0  13  ^^05 

Figure  5.  Forecast  of  the  surface  backscattering  image  from  the  MODIS  particle  “Seed”  Field  (hour  1).  The  24 
hour  forecast  of  the  particle  field  shows  the  change  in  MS  river  plume. 


4.  CONCLUSIONS 

Ocean  color  imagery  provides  a  realistic  observation  of  the  ocean  surface  optical  properties  but  only  to  a  depth  of 
approximately  1  optical  depth.  Optical  layers  many  times  extend  beyond  the  satellite  penetration  depth  and  are 
significant  different  than  those  observed  at  the  surface.  To  provide  a  more  complete  3-d  realization  of  the  ocean  bio- 
optical  properties,  we  have  demonstrated  a  method  to  couple  these  surface  satellite  observations  with  physical  properties 
from  ocean  models  to  extend  the  surface  optics.  The  procedure  is  based  on  a  characterizing  the  shape  of  a  Gaussian 
vertical  profile  based  on  the  interaction  of  the  mixed  layer  depth,  intensity  of  the  mixed  layer,  1  percent  light  level  and 
the  surface  satellite  chlorophyll.  The  method  used  observation  such  as  gliders  to  help  define  the  coefficients  which 
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describe  the  Gaussian.  The  example  of  the  3-d  profile  is  for  the  Northern  Gulf  of  Mexico,  however,  the  methods  has 
been  applied  to  other  areas.  The  methods  have  been  evaluated  with  insitu  observations  with  reasonable  results. 

The  forecasting  of  the  coastal  optical  properties  was  demonstrated  by  “seeding"  a  coastal  circulation  model  with  the 
surface  particle  distribution  and  concentration  based  on  MODIS  satellite  derived  backscattering  coefficients  using  the 
QAA  algorithm13,14.  The  intensity  of  the  backscattering  coefficient  (bb551)  was  used  to  determine  the  particle 
concentration.  This  particle  field  was  advected  on  hourly  time  steps  using  the  NCOM  circulation  model  using  a  fixed 
dispersion  and  particle  settling  term.  The  forecast  particle  field  was  accumulated  and  inverted  back  to  the  backscattering 
coefficient  and  compared  with  the  next  day  MODIS  backscattering  image.  The  difference  field  of  the  forecast  and 
satellite  image  provide  a  method  to  confirm  the  physical  circulation  field.  These  differences  additionally  illustrate  that 
simple  particle  advection  accounts  for  a  large  portion  of  the  coastal  bio-geo-optical  processes  occurring  with  24  hours. 
However,  the  difference  image  can  be  used  to  identify  other  bio-optical  processes  such  as  phytoplankton  growth,  particle 
upwelling  and  downwelling,  and  particle  fluxes  associated  with  bottom  resuspension  and  accumulation. 

We  provide  an  initial  effort  to  demonstrate  the  integrated  use  of  satellite  bio-optical  products,  numerical  circulation 
models  and  insitu  observations  for  monitoring  coastal  processes  and  provide  a  forecast  capability. 
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